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Abstract 



For first order differential equations of the form y' = J2 P =o Fpi x )y p an d 
second order homogeneous linear differential equations y" + a{x)y' + b{x)y = 
with locally integrable coefficients having asymptotic (possibly divergent) 
power series when \x\ — > oo on a ray arg(x) =const, under some further 
assumptions, it is shown that, on the given ray, there is a one-to-one cor- 
respondence between true solutions and (complete) formal solutions. The 
correspondence is based on asymptotic inequalities which are required to be 
uniform in x and optimal with respect to certain weights. 

1 Introduction and Main Results 

The main purpose of the present paper is to give, in terms of uniform asymptotic 
estimates, a precise meaning to complete asymptotic expansions (e.g., as power 
series followed by exponentially small terms) of solutions of a class of differential 
equations in a neighborhood of an irregular singular point (chosen to be infinity). 
The study of the exponentially small terms in asymptotic expansions has known a 
rapid development in the last years, especially after the pioneering works of Ecalle 
[TT] and Berry 
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For first order polynomially nonlinear (or linear) and second order homogeneous 
linear ordinary differential equations we give, using asymptotic inequalities, a one- 
to-one correspondence between formal solutions and true solutions. The represen- 
tation of a given function y turns out to be (under some assumptions) independent 
of the differential equation(s) of which y is a solution. 

The key ingredient is the concept of uniform optimal asymptotic inequalities 
which we illustrate in the following. 

Consider a formal power series / = Y^k=o fk%~ k , where fk are complex numbers 
and x is thought of being a large variable. We say that a function / is uniformly 
asymptotic to the series / with respect to the weight w, along a given ray in the 
complex plane, say 1Z = x > x > iff 



/(*) 



n-1 
k=0 



< w(n)x 



Wx G 1Z and Vn G IN 



Any function / that is asymptotic to the series / is uniformly asymptotic with 
respect to some w. The minimal such w is obviously 



w(f, /; n) = Wf(n) := sup x r ' 



X>XQ 



n-1 



f( X ) ~ J2 fk X ~ 



k=0 



1.2) 



(we might think of Wf as a transform of / with respect to /). With no restrictions 
on w we thus obtain the Poincare asymptoticity. With more conditions on w we 
obtain sharper asymptoticity classes: e.g. if fk ~ ^(nl) 13 and we consider the 
functions / for which Wf(k) < 7™(n!) /3 , 7 > a we obtain the familiar Gevrey or 
Gevrey-Roumieux classes [TS] . 

It is natural to take as a measure of the separation between a function / and 
a formal series, /. We will say that / is closer to / than g iff 



Wf{n) < w g (n) for all large n 
or, on occasions if a weaker condition holds: 



1.3) 



Wf(n) 

lim SUp ; ; < 1 



7WOO Wg{n) 



Given a formal series /, a ray 1Z and a class of functions J 7 , it is also natural 
to ask what are the sharpest bounds compatible with it, i.e. what is the "greatest 
lower bound" of the weights w such that w = Wf for some f. 
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We say that a function / is optimally asymptotic to / with respect to T along 
1Z, written 

oo 

/ ~ ^2 fkX~ k as x — > oo on 1Z (1.5) 

k=0 

and correspondingly that Wf is an optimal weight for / iff (J1.3)) (or sometimes (|1.4|) ) 
holds for a// g ^ T . 

As a first example, note that if /o is a convergent series and xo > > -R -1 
where i? is the radius of convergence, and T is the set of all functions defined for 
x > x then 

oo 

f-fo^f = fo = Y, hk X ~ k 

k=0 

Indeed, Wf (n) < Cp n (l — pxQ l \ — > as n — > oo and therefore uniform asymp- 

toticity of a function / with respect to w implies the convergence of the series f 
to /. " 

Next, take f\ = J2T=o k\x~ k ~ x . We will see that an optimal weight for this series, 
in the sense (jl.4j) . along the ray x > x > behaves like w(n) = (a*y/n-\- 0(l))n\ 
for large n, where a* is defined in (jl.l6|) . An optimally asymptotic function is 
f(x) = e~ x Ei(x) := e~ x P J^t^e'dt. 

It turns out that if fx, f% are two functions optimally asymptotic to f\ then 
fi — f2 = o(e~ x ) for large x. It is then meaningful to extend the notion of optimal 
asymptoticity to more general structures e.g., power series followed by exponentially 
small terms: in our example we can define 

00 fcl 00 A- 1 

/(*) - E -^TT + ^ /(^) - Ce- X - E -^TT (1-6) 

fc=0 x fc=0 x 

since the second relation can hold for at most one value of C . On the other hand 
e~ x Ei(x) + Ce~ x is the general solution of the equation /' + / = x~ l whereas 
/ + Ce~ x is the general formal solution of that equation. It follows that the relation 
f(x) — f + Ce~ x establishes a one-to -one correspondence between true and formal 
solutions for this equation. 

The main goal of the present paper is to show that by letting T to contain 
solutions of a class of first and second order differential equations, along some ray, 
optimal asymptoticity (in either the stronger form (jl.3j) or the weaker one ()1.4j) ) 
always gives a one-to one correspondence between true and formal solutions. 



3 



As a manifestation of the Stokes phenomenon, the complete expansion of a given 
solution of a differential equation will depend on the ray along which the asymptotic 
estimates are considered. In our example, we have 



e x £i(x) ~ 



3=0 



xi +l 



+ 




— ixi e 



—X 



X 



if arg(x) = 
if arg(x) G (0, 7r) 
if arg(a;) G (— ir, 0) 



(1.7) 



For a complete study of the Stokes one would need to consider besides rays, 
parabolas (of the form arg(x) = A | a: | ~ x / 2 for fll.7|) ); then the constant beyond all 
orders changes smoothly in a narrow region near the Stokes line (it is a smooth 
function of A in the example above), see |14j . 

However in this paper we are not making analyticity assumptions on the coef- 
ficients of the differential equations and we will be consequently mainly concerned 
with the behavior of solutions along a fixed ray. 

A function is characterized by its generalized expansion with a precision compa- 
rable to that of truncation to the least term of the complete asymptotic expansion. 
The relation with this technique is explored in section 2. 

In this respect we mention the paper jTH] in which it is shown that truncation 
to the least term can be used to measure the terms beyond all orders for second 
order linear homogeneous differential equations with coefficients analytic at infinity. 
The corresponding results in our paper show that the same is true for equations 
with not necessarily analytic coefficients (and as such only defined on a ray) having 
possibly divergent series at infinity (divergent no faster than factorially) and first 
order polynomially nonlinear equations. Removing the analyticity and linearity 
assumptions raises difficulties requiring new techinques for the proofs. In particular 
we provide a method of estimating the growth of the coefficients of formal series 
solutions based on the recurrence relation that they satisfy. We mention in addition 
the paper ^7] in which very interesting re-expansions for the optimal remainder are 
obtained in the analytic case. 

We also discuss the possibility of extending our results to more general systems 
of differential equations. 

We postpone further discussions and examples until after Theorem 11.11 below. 
In Section 2 we study, in some generality, the connection between optimal asymp- 
toticity and the method of truncation near the least term. Section 3 is devoted 
to the proofs and in Section 4 we discuss, in the context of a particular family of 
differential equations, the question of typicality of the divergence of the asymp- 
totic representations and give a decomposition formula which provides a different 
perspective on complete asymptotic expansions. 
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Consider the second order linear differential equation 

y" + a{x)y' + b{x)y = (1.8) 

given on the ray 

TZq = |x : e~ td x > x \ 

The coefficients a(x) and b(x) are assumed to be in L} oc (7Z) nd to have asymptotic 
power series at infinity of the form 

oo oo t 

flW~ES ; 6 (*)~E* (^-+oc) (1.9) 
fc=0 X fc=0 x 

Moreover, we require that the functions a(x) and b(x) satisfy a Gevrey-like 
condition JH] on the given ray namely, for some k < \oq — 46 | _5 



a n , b n < constK™n! 

\x\ 



n— 1 _ 

l a ( x ) — E — fel < C 0118 ^™^! 



fc=o x 

n-1 



IK^-E^ll < constK n n\\x\- n (1.10) 



fc=0 

uniformly in x G TZg ,n E IN. 

If the conditions (jl.K)j) hold, then the rate of divergence of a formal series solu- 
tion depends only on the first few terms in the asymptotic series of a(x) and b(x). 
The results below (not merely the proofs) depend on this assumption. 

The polynomial 

A 2 + a A + & (l.H) 

is assumed to have distinct roots, A 1; A 2 (if the roots coincide then there are no terms 
beyond all orders to worry about). The following expression is a formal solution: 

S = de XlX x ri S^x) + C 2 e X2X x r2 S 2 (x) (1.12) 

where for 2 = 1,2 

d\ \ + b\ 
a + 2\i 
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and 



fc=0 x 



are formal power series. The asymptotic behavior for large n of the coefficients of 
the power series, as follows from the Proposition 13.21 below, has the form 

^-r+g,— J (a, -a,)- (L13) 

(the expression for s 2)n is obtained from the one above by interchanging the indices 
1 and 2). Choosing 

Si ; o = s 2; o = 1 

the coefficients of the series are uniquely determined. 

The expression (11.12)) containing two arbitrary constants is the general formal 
solution of our equation, in the differential algebra generated by power series and 
exponentials (see 0, |Zj, [HP- 

The notion of optimal asymptoticity can be extended in a natural way to asymp- 
totic structures of the form ()1.12|) . We say that / is uniformly asymptotic to S for 
a weight (w 1 ,w 2 ) along the ray Tig iff 

y(x) - d x ri e Al x V ^ - C 2 x r2 e Aa x V ^ 



< (1-14) 

|s ri e AlX |iwi(A;i) + |x r2 e A2a; |w 2 (fc 2 ) Va; G 7l e and fci,A; 2 G W 



We say that y is optimally asymptotic to S along the ray Tig and we write 



y~Cix ri e Aia; y-^ + C 2 x r2 e A2:c y ^ (1.15) 

if there is a weight (wi(y; n), w 2 (y; n)) such that f)1.14j) holds and moreover for any 
/ G Wi(y; n) < w\(f\ n) and w 2 (y; n) < w 2 (f; n) for all large n. 

Theorem II . li b elow states that for each formal solution there is, on the given ray, 
a unique true solution of the equation which is optimally asymptotic to it relative to 
the space of solutions of the same differential equation. The theorem can be easily 
extended to encompass differential equations that can be brought to the form that 
is treated here by an algebraic change of independent variable and linear changes 
of the dependent variable; see also Remark 4 below. 
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Let B be a (large enough, positive) constant, r := r 2 — r 1; p = 3?(r), and 
A = A 2 - Ai. 



Theorem 1.1 For any Ci, C 2 there is a unique true solution y(x) of the differential 
equation M.ty) . under the assumptions following it, with the property \1.14\) for the 
weight 



with 



Wi (y; h) = IdKAJki + B)T(h - (-l)»|A|- fe ' i = 1,2 



a\ Rl \ ifetfr-Ve* 
(J otherwise 



where 

a* : 



V^maxe-* 2 f e t2 dt = 0.765151.. (1.17) 

x>0 Jo 

and A 2 is obtained from 11.16)) by interchanging the indices 1 and 2. 



Comments. 

1. ) Let y(x) be the unique solution provided by the theorem and choose the spe- 
cial truncation orders k± = k 2 = [\xA\], which correspond generically to "truncation 
to the least term" of the series. It will follow from the proof that the remainder (lhs 
of (j!.14|) ) is in this case less than 

const -^ e ~ lxAl max { \x r2 e XlX | , \x r ^e Mx | } (1.18) 
J\x\ 

i.e., asymptotically smaller than any (nonzero) solution of the equation by a factor 
of J\x\ on the Stokes line S(xA) = and by an exponential factor in the generic 
case $s(xA) 7^ 0, indicating why there is uniqueness of the representation. There is, 
in effect, a connection between the shape of the optimal weight and the precision 
of the least term truncation method, as will become clear in the folowing section. 
Moreover, the solution that is best approximated by the simultaneous least term 
truncation of the two component series in ()1.15|) is precisely the one having the 
optimal weight. There are differences between the two approaches but the terms 
beyond all orders that they predict agree for this class of differential equations. 

2. ) The values of the constants A\ and A 2 in ii) are crucial for the result. If 
say, ^(e' 6 ^) < 3?(e l6l Ai) and Ei / then there will be no solutions satisfying 
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the required inequalities with Ai < a*\Ri\ whereas if A\ > aJRi\ there will be 
infinitely many such solutions. In contrast, the theorem is true for any large enough 
B. We also want to stress that uniqueness is relative to the space of solutions of 
the differential equations treated here. 

3. ) As a manifestation of the Stokes phenomenon ([2], [E| 5 the constants 
Cj in the asymptotic representation of a given solution depends on the direction 9 
of the ray considered (provided, of course, that the differential equation satisfies 
our hypothesis on more than one ray at infinity). The example f)1.7|) given at the 
begining of the section illustrates this point. 

The derivation of ()1.7|) is given at the end of Section 2. We also mention at this 
point that a more detailed analysis is possible along the same lines and it shows that 
the term beyond all orders varies smoothly on a scale of the order arg(x) ~ x -1 / 2 in 
a way which agrees with the results obtained with the hyperasymptotic technique 
of Berry [Tlj, We will however not pursue this issue here. 

4. ) One can in principle allow for more general formal structures than (jl.l2j) . for 
instance those obtained by substitutions and formal algebraic operations on (jl.l2|) 
(for a discussion on what formal structures are relevant to solving differential equa- 
tions see [ini), to allow for asymptotic representations of the solutions of equations 
that are not of the form required by the theorem, but can be brought to that form; 
the asymptotic inequalities are well suited for simple algebraic operations. Consider 
for instance homogeneous Airy equation 

y "-xy = (1.19) 
After the substitution y(x) = exp(| x^)g(x) followed by taking x = s 2//3 we get 

9 " + ( 5 + S )9 ' + S 9 = (1 ' 20 > 

to which Theorem 11.11 applies and we obtain, after undoing the transformations, 
the following asymptotic representation of the general solution of the Airy equation 

3 00 3 00 / \k a 

y(x) * C x eh- X -V± £ ^ + C 2 e~^ £ (1.21) 

k=0 x 2 k=0 x 2 

when x becomes large along a given ray in the complex plane. For a fixed solution 
Ci, C2 will depend on the ray. A brief derivation of (jl.21j) and the expressions of 
Ai, A 2 and Sk are given at the end of Section 2. 
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5. ) Whenever e l0 A ^ IR _ all the terms of the expansion (J1.12)) are simultane- 
ously visible to the inequalities, as it will follow from the proof of the theorem. If 
on the other hand e l9 A < i.e., on the Stokes line, besides the first series only the 
first term of the second series is caught by the inequalities. However this first term 
is enough in order to establish a one-to-one (linear) correspondence between formal 
and true solutions, so that we keep, by courtesy, all the other terms of the second 
series. 

6. ) For obtaining complete asymptotic representations by means of optimal 
inequalities for higher order differential equations or systems of equations it seems 
that it is necessary to use, inductively, asymptotic comparisons with solutions of 
lower order ODE's. 

The result below shows that a given function cannot be a solution of two es- 
sentially different second order differential equations that satisfy our assumptions. 
We say that a differential equation of the type considered in the theorem is in the 
canonical form if: 

6 = 0; Ai = 0; |oo| = 1; n = 0; b = h = (1.22) 

This can be always achieved by the substitution y(x) = exp(Ai x)x ri y(x) followed 
by a change of independent variable x = e ie \\i — \2\~ 1 x. The differential equation 
is said to be formal if the coefficients a(x), b(x) are replaced by their formal series. 

The condition (|l.l(Jj) plays a crucial role in the result below. 

Proposition 1.1 Let S = J2^Lo s nX~ n be a formal power series such that s n have 
the asymptotic behavior M.l'Jji with R\ ^ 0. Then S formally solves at most one 
second-order canonical formal differential equation of the form if the formal 

series of a(x) and b(x) satisfy the first condition in il.l(J\) . 

It follows for instance that within the class of differential equations we are con- 
cerned with, there is a one-to-one correspondence between generic formal solutions 
and true solutions. Indeed, a true solution is associated with a unique representation 
within the solutions of the differential equation it originates in. But then, the rep- 
resentation determines uniquely the power series of the coefficients of the equation 
and by this determines uniquely the differential equation itself. Conversely, a for- 
mal expansion solution is associated with a unique differential equation and within 
that differential equation with a unique true solution. In this sense we also obtain 
a summation method. One can attempt to continue this construction inductively, 
considering at the n-th step differential equations whose coefficients are solutions 
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of the equations gotten at the step n — 1 and obtain a more general representation 
spaces and summation methods. 

* 

Similar results hold for first order nonlinear differential equations of the form 

y' = F{y,x):=Y,F p {x)y* (1.23) 

p=0 

defined on the ray IZg. 

We require that for p = 0, .., P, F p (-) G C[R e ], F(0, oo) = / , = 0, F y (0, oo) =: 
a = — /x,o 7^ and impose the Gevrey-like condition 



n—l r 

Jp,k 

k=0 



If, 



p,n 



< C F K n n\\x\- n 

< C F K n n\ 



;i.24) 



(where Cp is a constant) uniformly in x G TZg, n G IN, where k < 

We are only interested in the case when the equation allows for exponentially 
small terms beyond all orders; the condition for that (as it will become immediately 
clear) is 



arg(xa) G (--, -) 
Lemma 1.1 There exists a formal series solution 



k=l 



X h 



;i.25) 



;i.26) 



of the equation The coefficients of the series have the asymptotic behavior 

(oo p \ 
R + E -^U~ kT ( k + r ) (zsk^oo) (1.27) 

where 



r = fi,i + 2/2,0/0,1/tt 
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If we are now looking for nearby formal solutions Sg := S + S we get 



6' = l-a+ 1 -^)5 + 0(5 2 ,5/x 2 ) (1.28) 



x 



so that for large \x\ 



5~Cx r e~ ax (1.29) 

The analog of Theorem 1 1.11 in this case reads: consider the equation (|1.23|) under 
the assumptions ()1.23|) - (jl.25|) and let p = 9R(r). Let B be large enough. 

There is a one-to-one asymptotic correspondence between true solutions which 
decay at infinity and formal solutions: 

Theorem 1.2 Given any constant C there exists a unique true solution y(x) of the 
differential equation such that 



k-i 

y{x)- — ~Cx r e~ ax 



for all x G IZg and k G IN, where 



< (AVk + B)T(k + p)\\x\~ k (1.30) 



A= [ a*\R\ ife i9 a G IR + ^ 
1 otherwise 

Conversely, given a solution of such that xy(x) is bounded for large x, 

there is a (unique) constant C such that, for all x <E7Zg and k G IN, il.3U\) holds. 

In the sense of (jl.30|) we then write 

oo e 

y(x) ~ Y. — + CxTe ~ aX ( L32 ) 

3=0 Xl 

The comments that follow Theorem ll.il with obvious adaptations, also apply in 
this case. 

We could actually continue the construction of the formal solution and consider 
the complete formal solutions or "transseries" ( >C3) which in this case have the 
form 

Y + x r e- ax Y x + x r2 e~ 2ax Y 2 + ... 
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in which Yj are formal power series. Y\ is determined up to an arbitrary constant 
which, once given, determines completely all the following power series Yj, i > 2. 
We will however not pursue this direction here but merely remark that the one-to- 
one correspondence between formal solutions and true solutions is pinned down by 
the first term of the second formal series and contend to control the asymptotics 
only to that level. 

2 Optimal estimates for a class of divergent series 

In this section we estimate the optimal bounds for series that diverge in a way 
typical for differential equations and make the connection between optimal uniform 
inequalities and the method of truncation to the least term. We also give results on 
the precision with which a function is represented by its optimal asymptotic expan- 
sion. We consider formal power series, for large argument on a ray and establish 
the connection to the technique of optimal truncation of divergent series. 

fc=0 x 

where 

Cfc = (1 + e k ) k i4) C k exp(G(Jfe)); e k -> as k -> oo (2.2) 

Since the term C k can be always absorbed into the independent variable so that we 
will take C — 1. 

In (J2.2)) is real and we make the following assumptions. G(-) G C 3 ( 1R + ) 
is a real- valued increasing, convex function with the properties G"(x) — > and 
x 2 G"'(x) — > —7 < as x — > 00. Such a series is necessarily divergent for all values 
of x (convergent series are discussed in the introduction). An example satifying the 
requirements would be c k = C k Y^{k + r). 

The behavior of the optimal weight for a divergent series depends critically on 
the ray along which we consider the asymptotic series. In the setting (|2.2j) the 
optimal weight is larger by a factor of order ^/n along the real positive axis than 
it is along any other ray. This is a manifestation of the Stokes phenomenon at 
the level of the asymptotic series themselves. Intuitively we can account for this 
behavior in the following way. A natural scale for studying the difference between 
a function and the n — th truncate of its asymptotic series is the n + 1-th term 
of the series (indeed, the n + 1-th term is meant to be a correction for the above 
mentioned difference). When x is very large, the successive terms of the series start 
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by decreasing fast as so does the error in approximating the function by its series. 
The least error suggested by our rough guide, the next term of the series, reaches a 
minimum when 



The width of the minimum is of the order y/s^ and within this width, for x > 0, 
the ratio of two terms is approximately one. But then we realize that if at one point 
within this range the difference between the function and the truncate of the series 
is approximately equal to the next term there will be a point within the same region 
where the difference will be roughly y/n times as big. If x is not real the ratio of 
the successive terms is of the form e~ %e so that the overall accumulation of errors is 
still of the order of a constant. 

We want to stress that, as it follows from Lemma [3.11 below and the estimates 
in its proof, that the least-term truncation of a formal series solution is a good 
approximation only for one solution of a given differential equation; for all all the 
other solutions there are exponential corrections that are much larger than the least 
term of the series and have to be removed before calculating the function from its 
series (to such an accuracy). 

Consider the power series S, under the given assumptions, for x on the ray TZg. 

Lemma 2.1 i) For any 9 there exist (smooth) functions that are optimally asymp- 
totic to the series \2.2\) . A optimal weight is 



More precisely, assume that $ is a function with the property that 3r] and B such 
that V x G IZg and Vn G IN we have 



G'(s x ) = \n(\x\) 



(2.3) 





(2.4) 



Then, r] > 1 and there exist functions satisfying \2.J\I for rj — 1, for some B. 
ii) For 6 ytz assume further that 



(2.5) 



Then the optimal weights can be taken 



g*(ri) = (a(6) + Bn )\c n+1 
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where the constant a{9) can be computed explicitly and for small 9 has the behavior 

a{6) ~ 9~ l (2.6) 

The behavior \2. 6]) is actually present in the concrete estimates of the exponential 
integral, see \3. 83\) . 



The result below answers the existence part of Lemma \2. 11 i) and gives the connec- 
tion with the technique of summation to the least term. Let n x be defined, for x 
large enough by n x = [s x ], (cf. (|2.3|) : [ ] denotes the integer part) 

Proposition 2.1 i) Let 9 = 0. A function $ satisfies the uniform inequalities \2.Jj!{) 
with r\ — 1 iff 



lim sup n. 



-1/2 



X 



$(x) 



n x ,, 
rJ 







(2.7) 



Instead, if the limit above is e > and $ satisfies \2.4\j for some r\ then rj > 1 + e. 

ii) If 9 7^ and \2. 5)) holds, then a function $ is uniformly asymptotic to the 
series \2.2}) along the ray TZg with respect to a weight 

win) = C\c n+l \ 

for some C iff 



lim sup 

\x\^oo 



X" 



j=0 



X 3 



< oo 



(2.8) 



Corollary. If $ 1; $ 2 are two functions satisfying the hypothesis of part i) of 
Lemma f2. II with rj = 1 then 

limsupn^ 1/2 r r|$i(ar) - = 

x-i-oo I C nx | 

In particular when G(n) = ln(r(n + n )) then 

$i(x) - $ 2 (x) = o(e~ x x no ) (2.9) 
Correspondingly, in the case iii) we have 
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limsup^ 1 / 2 -^ r|<&i(aO - ®2(%)\ = 

x^oo \ c n x \ 

and for G(n) = T(n + n ), 

$ x (x) - $ 2 (x) = (V N x n °- 1/2 ) (2.10) 

Classically, two functions $i and $2 have the same asymptotic series (provided 
they have one) if their difference is asymptotically less than any power of x. The 
corolarry shows that there is a definite precision gain when requiring optimal uni- 
form estimates. 



3 Proofs and further results 



Before giving the general proof we mention as an illustration the particularly easy 
computation of the optimal weight for the exponential integral: see the arguments 
starting with eq. (j3.8(ij) . 

In this section we make the following convention: for n > k, 



E 



E 



"3 ~ 
j=n j=k 



a ; 



For the proof of Lemma f2. II we need the following elementary result. The nota- 
tions and hypothesis are those that precede Proposition 12.11 



Proposition 3.1 i)Let ko be large enough. Then for 6 = 0, 

K 



lim sup I (K + k )- 1/2 x K e- G{K) 
™ KtIN( 

The supremum in \3. 1)) is actually attained for 



E 

j=n x 



El 

X J 



a*7 



(3.1) 



K ~ n ±:= [s x ± 



1 [Y 



11) J/0 ^ then 



lim sup < x K e~ G[K) 
™ KeM { 



K 

E 

j=n x 



El 

X J 



< 00 



(3.2) 
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Proof 

The proof is essentially straightforward. Note first that from the assumptions it 
follows that G"(x) ~ j/x for large x. 

Note also that we need only consider the case 6j = 0. Indeed, the presence of a 
(large enough) constant ko makes the proposition above insensitive to the behavior 
of Cj for small j. On the other hand, 



K c- 
A xi 

3=n x 



K 



i4> p ~G(j)+j ln(x) 



3=n x 



K 



+ 0( max | ei |)) J2 e 



-G(j)+j ]m(\x\) 



K<j<n 



(3.3) 



3=Tlx 



and the maximum above approaches zero as K — > oo. 

• i) Take first = 0. Let (3 G (|, |) and start with the range of .fT so that 

\K -n x \ < nf 

which is actually the important range with respect to (|3.1j) 

Here, the Euler-Maclaurin summation method is suited. For definiteness we take 
K >n x (the other case is very similar). We have 



^_ e G{j)-j\nx _ + Q^ s W-2^ e \G"{s x )(j-s x f 
C n x 

(l + 0(s 3 /- 2 )) [ +l eho"{s x )U-s x ? dt = 
(l + 0(sf~ 2 )) [ j+1 e&'^W-^dt 



(3.4) 



so that 



max 

\K-n x \<n% 



K -l/2 e -G(K)+K In x e G(j)-j In s 



(l + 0(nf - 2 )) max J jr-^e"^^--)' T +1 e^"M(*-^) 2 rft I = 

V ' \K-n x \<n% { Jn x J 

(l + 0(nf- 2 )) a*(KG"(s x ))- l i - (3.5) 



as x — > oo. 
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• It remains to obtain an upper bound in the region \K — n x \ > nf. Assume 
K < n x (the opposite case is treated similarly) and let q = s x — |s^. With F(z) : = 
G(z) — z ln(x) we have in view of the estimates above, 



e FU)-F(K) < e m-HK) + co ^t^de F ^~ F ^ (3.6) 

j=K j=K 

On the other hand, from the definition of s x , 

F(j)-F(K) = ^ dz f dtG"(t) < const f* \n(z/s x )dz 

JK Js x JK 

Throughout the domain of integration s x > z + > z + z 13 so that 

ln( , /Si)<ln ( 1 __i!_) < _I 2 ^ 



and we get, for some positive constants C±, C2, 

j- e F(j)-F { K) < ^ e -C l[{ K + n)^K^ < ^(Jf + ^l-fl < \ W±S^ 
K n=0 ^ V 7 

In the same way, for some positive constants, 



e F(q)-F(K) < e -Ci[qf>-K?] < ^ChK^q-K) < g-Csg' 3 < e -C&Ph 

(3.7) 

so that the second term in (|3.fi|) vanishes as x — > 00. 

• If now 7^ the upper bounds are trivially obtained by taking absolute values 
in the sums. For a lower bound we note simply that 

n+ „ K 

E £i = n? E e -^i K.) (l+0(n^ 2 )) 

• ii) Let again = G(A;) — A;ln(|x|) — z"&0. Given a function /(x) which is 
asymptotic to the series S we write 

m = : E^- i^xo*) ( 3 - 8 ) 

A;=0 X 
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and we will choose x to make f(x) optimally asymptotic to S. Consider first the 
region \k — n x \ < and define 



fcirt> e F(k) 

'e~ ie - V 



7e 



i(k—n x ) K—n x —\ 



-tj$ 





3=-nz 



It is easy to check that 



Cfa / , 1/9 

(Tk+l ~ cr k = — (1 + o{n x 1 



(3.9) 



(3.10) 



(to get ()3.9j) . it is simpler in this case to solve ()3.1()|) directly by perturbation 
expansion than to use the Euler-Maclaurin summation formula). Using ()3.8jl . ()3.9j) 
and ()3.10|) we get 



, -i6\ 



J(n x -n)6 



X 



n+1 



c n+l 
i9h x (n) 



k=0 



Ck_ 

X k 



X + o(l))e~ 7 2 ™z 



(3.11) 



where ^(n) = 1 if n > n x and is zero otherwise. For large x ()3.11|) can be visualised 
as distance between z = 1 to the points on two spirals centered at e 



-idh x (n) 



X: 



and with diameter slowly decresaing as \n — n x \ increases. It is not hard to see that 
there is in this region a best choice x* °f xi x ) (which minimizes the maximal error 
in ()3.11|) since for any value of a; there is only a finite set of n such that |n— n x \ < . 
For 6 incommensurate with tc the geometry of the problem shows that 



X* = e 2cos- + o(l) 



(3.12) 



for which we get 



a(9) 



1 + | sin || 



sin 



(3.13) 



If 9 is a rational multiple of it the computation can still be done explicitly, in a 
straightforward manner and a{9) is slightly less than (|3.13j) . 

In the region | k — n x \ > a similar calculation shows that the error is bounded 

by 

E n <\l- e *W|-i + (l)| x ( x )| < a{9) + o(l)\ X (x)\ 
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which shows that x* is indeed the optimal global choice of %, as it is not possible to 
decrease the error in this region without making it unbounded in the n — n x \ < rzf 
region. 
□ 

Proof of Lemma 12.11 For part i), if we assume there existed a function <£> sat- 
isfying the estimate (|2.4j) uniformly in x and n for some r) < 1 it would follow 
that 

±(^)-E|) <Va^-HV^ + B) C ^ (3.14) 

By adding the two inequalities above we arrive at an immediate contradiction with 
Proposition 13.11 

Part ii) and iii) are straightforward consequences of Proposition 12.11 Indeed, 
note first an easy example of a function satisfying (|2.7j) and ()2.8j1 : 

*•(*):= £i (3.15) 

3=0 X 

It is not dificult to smooth and preserve (|2.7|) . (|2.8|) since the size of a jump 
at a point of discontinuity of is |c n:E x~ na: |. 

On the other hand, Proposition 12. ll follows immediately from Proposition 13. ll bv 
the triangle inequality. □ 

* 

Proof of Theorem 11.11 

There are no assumptions that would distinguish between the quantities with 
subscript "1" from those with subscript "2". The existence part of the theorem 
follows (trivially, using triangle inequalities) if we prove it for C\ = 1 and = 0. 
This case will follow from Lemma 13. II below. 

We can assume that the differential equation is brought to its canonical form 
because the asymptotic inequalities are transformed in an obvious way in the sub- 
stituions involved. Also we note the following inequality that we will use frequently 
and which is clear from the integral representation of the Gamma function: 

|r(z)| < r(»(x)) for ft(x) > -1 

Let r := r 2 = — a\, p := 3?(r) and a := oq. 
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Lemma 3.1 There is a solution Y of the differential equation U.8\) under the as- 
sumptions il.y\) - il.2ty) which satisfies the inequalities: 

n 

\Y - -^r\ < const e~ x x p ~ 1/2 for xe(n,n + l) (3.16) 

fc=0 x 

where const does not depend on n, x 

Before we prove this Lemma we need some results on the formal solutions. Let 
Sk = si f fc. The condition that J2kLo fr is a formal solution of the equation leads to 
the following recurrence relation for the coefficients of the series: 



s n = - [n - 1 + r + — — J s n _i + V — (U - n)aj + b j+1 ) s n _j (3.17) 
Proposition 3.2 The behavior of s n for large n is 

s n ~(R+f:^)a- n r(n + r) (3.18) 

\ m=l 71 ) 

In the assumption of analyticity at infinity of the coefficients of the equation, 
the leading behavior of s n was shown in [TO] . 

Remark. Although in a typical case the constant R will be nonzero as it is easy 
to understand by examining the recurrence (|3.20j) below, massive cancellations are 
possible so that the series R+^ + .. could be zero to all orders in 1/n; this case is 
important in its own right-see Section 3 for some examples and applications in this 
connection. 

Proof 

In order to avoid dealing with the poles of the T function which might occur for 
(uninteresting) small values of n we let 

r r r(x) ifK(*)>-i (319) 

x ' | 1 otherwise 
It is convenient to pull out the leading behavior suggested by ()3.17|) : 



s n . ct r ra _i_ r c n 



We get for c n 
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c„ = c n _! + +0(n- 2 )c n _ 2 + ^ C(n,i)c n _ i + o(n- 2 ) (3.20) 

where the term with j = 3 and n > max{p, 3} terms with j near n in the sum 
were treated separately. Due to (jl.lOjl the latter give a collective a contribution of 
o{n~ 2 ) and also, the following estimate holds for C(n,j): 

\C(n,j)\ < const ^—IM = 0(n- 3 ) (3.21) 
n! 

in the range n — 3 > j > 3 (the inverse of the binomial coefficient is convex in j). 

By induction on n we see that that \c n \ < A Y[j <n {\ + const j~ 2 ) for some 
large enough A and const, whence the sequence c n is bounded. But then it follows 
immediately from ()3.20|) that the sequence {c n } is convergent. Furthermore, taking 
c n = d n Y[j< n (i + Aj~ 2 ) we see that 

d n = c?„-i + -J (3.22) 
where is a bounded sequence and thus, if R is the limit of the c n we have 

c n = R + Oin- 1 ) (3.23) 

It is now easy to bootstrap the estimates for c n in the recurrence (taking out ex- 
plicitely more and more terms from the sum) to get c„ = R + — + 0{n~ 2 ) and so 
on. 
□ 

Let 

U n ={\R\ + 0(n- 1 ))Y(n + p) 

represent an upper bound of \s n \. Note that if R ^ then \s n \ ~ U n for large n. 
Let 

fc=0 X 

Take n to be large enough so that for x > n the asymptotic estimates for a(x), b(x) 
hold and so that n + p > 0. 
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Proposition 3.3 In the differential equation verified by Y r 

Y£ + a(x)Y^ + b(x)Y n = R n 
the inhomogeneous term satisfies the estimate 



(3.24) 



R n (x) = (R+Oin- 1 ))^"- 2 1 J n+2 ' ~ {V2^R+0{x- l ))a- n - 2 x r -^e- x (3.25) 
for x G (n, n + 1). 
Proof 

We substitute a(x) = X)o ffc + (and the corresponding expression for b(x)) in 
the expansion of the LHS of 1)3.24)1 . To simplify the notations we let throughout 
this proof a n+ i = A n+ %(x); b n+1 = B n+ i(x). The coefficients of x~i with j < n + 1 
vanish by the definition of Y n . The coefficients of x~ n ~ x add up to a(n + l)s n+ %. 
Indeed, — a(n + l)s n +i is the only term of order n + 1 which is missing in (J3.25)) 
with respect to the corresponding expression of Ym, M > n. So, 

Rn = a {n+ H S 2 n+1 + £ x-T m (3.26) 

X m>n+3 

with 

T m := d k - 1 ) s k-iO-j + s k bj) (3.27) 

k+j=m 
k,j<n 

An individual term in T m can be bounded by: 



| (A; — l)sk-i(ij + Sk bj\ < const U n 
as it is easy to see using estimates for Sk, a^, bk- Thus 

x~ n - 3 T n+3 < const n~ n - l U n 

and the same is true for J2 m > n +3 x ~ m T m , as shown by a gross majorization by the 
number of terms (0(n) by (|3.27)l ) times the maximal term. It follows that the 
modulus of the sum on the RHS of ()3.26|) is at most 

0(n->) T{n + r + 1) 
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for x G (n, n + 1). The rightmost estimate in (J3.25|) is just Stirling's formula. □ 
We return to Lemma 13. 11 
Proof 

The equivalent vectorial equation is now preferable: 

F' = AF (3.28) 

where 

•<: i ." ' i: T:-< ! 



-b(x) -a(x)J ' \f 
On the interval [n, n + 1] we look for solutions in the form 

F=(^)+H(rv) 

H has to satisfy the differential equation 

H' = M + T n ; T n = ^ _^ j (3.29) 

and the continuity condition (which is actually enough to ensure the smoothness of 
Y) 

H(n;n+) = H(n - l;n_) - E n (n±:=n±0) (3.30) 



-1 



A differential equation of the form (jl.Hj) under the given assumptions admits 
always two special solutions which for large \x\ have the behavior (see [S]) 

yx(x) ~ 1 + — + 0(x- 2 ); y 2 (x) ~ a; r exp(-a s)(l + O^ 1 )) (x -> oo) (3.31) 

It is convenient to choose a particular fundamental matrix of the system ()3.28|) 
constructed with these two solutions: 



MM .-{Vi(*) V2(x)\ fl + Oix- 1 ) (l + 0(x- l ))x r e-«* . 

W " l .^(x) 2/2^)/ V 0{x- 2 ) -a(l + 0(x- r ))x r e- ax > ; > ' > - ; 



and write the solution of ()3.29j) in the form 
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U(n,n + z) = M(n + z) QT M'^n + t) T n (t)dt+ M _1 (n) (H(n - 1, n_) - E n )) 

(3.33) 

Taking z = 1 in the relation above we get a recurrence relation for the {H(A; — 
1, k_)} k . With the substitution M~ l (k + l)H(fc, (k + 1)_) = q fc it reads 

qn = q„_i + / M-\n + 1) T n {t)dt - M _1 (n)E„ = q n _! + M" 1 (n) V(n) (3.34) 

JO 

with 

V(n) := / M _1 (n;7i + t)T n (t)dt-E n (3.35) 
Jo 

and where M(n; n + t) is the fundamental matrix specified by M(n, n) = J, whence 

H(fc,(fc + l)_)=M(fc + l) f J] Af-^VOO + qfco) ( 3 - 36 ) 
For large n and t < 1 we obtain in a straightforward manner, 

M- 1 (n,n + t) = (J t) )+0(n- 1 ) (3.37) 

and from (|3.25|) 

T„(n + t) = (1 + Oin- 1 )) ( ) i?» (3.38) 

Beginning with this point, most estimates will be different in the special case a = 1 
which corresponds to being on the Stokes line. Combining ()3.38|) and Propositions 
E21 Owe get from (j3~55j) 

V(n) = it^(n) ((^ ) + ( - 1 ) + OCn" 1 )) = ) i^n) (a = 1) 

(3.39) 

and similarly 

V(n) = 0(l)R n (n) («^1) (3.40) 
In view of (13.32)1 a direct calculation gives for M~ 1 (n) 
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\ + 0{n- 1 ) a^ + Oin- 1 ) 

,0(n~ 2 )n- r e an (-o^ 1 + 0{n- l ))n 

(note that the absence of exponential factors from the first row is not an effect of 
an approximation; in fact arg(a) is arbitrary and nothing is assumed about the size 
of the exponentials). Thus 

0(1) \ „ , v / Oirf-^e 



M-\n)V(n) = ( 0(n V e n } ) Rn(n) = {0^1/2^-1^ ) (« * 1) 



and correspondingly 



HV(n)=( vConstn _ r _ le „ji?» = ( v 0(n _ 3/2) j (a = l, 



so that the series Xlfc^ 1 (/c)V(fc) converges and 

0(n r - 1/2 i 

G ( n -l/2 e («-l)n) 



E^ 1 («)V(n) = (^i"^ ( tr)L) (3-41) 



k=n 

and 



If we make the choice 



in ()3.36|) we get 



( n -i/2) ) (" = ■'■) (3-^) 



q fc0 = - J] M-^fcjVC*) 

fc=fc + l 



H(fc, (k + 1)_) = —M(k + M-\j)\{j) (3.43) 

fc+i 

so that in view of (f!OTjl . (j3~4^jl and (jO^jl 

H(n, {n + 1)_) = 0(n r - 1/2 e- n ) (3.44) 

for all a with |a| = 1 □ 

At this point we can prove the existence part of Theorem 11.11 For R 7^ it is a 
direct consequence of the Lemma f3. 11 and of the Propositions 12 . II and 13. 21 If R = 
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we proceed in essentially the same way. Taking n such that x G [n, n + 1] we get 
for the special solution provided by Lemma [3. 11 



\Y 



3=0 X 



Y, 



+ \T— I 

j=fc x 



(3.45) 



Since now we have \s n \ < constr(n+p— 1)| the second sum in ()3.45|) can be estimated 
using Proposition 13. II taking Cj := constr(n + p — 1) (and G(j) := ln(c.,-)). It follows 
immediately from (|3.45|) that Y satisfies the inequalities (|1.14|) with G\ = 1, C 2 = 0, 
A\ = and some B. 

Uniqueness. We can assume without loss of generality (jl.22j) and 3?(« < 0). The 
case $l(a) = is actually trivial since it reduces to a statement about asymptotics 
to all orders. We show that if y satisfies (I1.14J1 then y = C\y\ + Ciyi where y\^ 
are the solutions corresponding to (Ci, C2) = (1,0) and (0,1) respectively whose 
existence has already been proven. 

We have to distinguish the case A\ 7^ (which implies a = 1 and R 7^ 0, 
cf. (11.16(0 . Taking k 2 = (which means that we only keep the leading term of the 
second power series) it follows from (|1.14|) that 



|tfi|(,Wtfi 



F(x)-C 1 J2^f-C 2 x r e 



B) 



x 



-k-X 



3=9 



X 3 



r(A; + r + l) + const \C 2 



< 



X 1 " 1 e~ x 



The case C\ — is trivial because the only solutions that decay exponentially 
at infinity are multiples of y 2 - With C% 7^ an elementary computation shows that 
one can choose a B' so that for all x G 1Z$ and k G IN 



dKAiVk + B) \x- k ^T(k + r + 1)1 + const \C 2 
\d\iAxVk + B') x- k - x Y{k + r + 1] 



x r - x e~ x 



< 



Consequently, the function y(x) — C 2 x r e~ x is optimally asymptotic to the series 
C\S\ (since 0). Because Ci?/i(a;) has the same property, we have, in view of 
the corolarry to Proposition 12.11 



\y(x) — C 2 x r e x — Ci yi(x)\ = o(x r e x ) = o(y 2 ) as x 



00 



which means that y(x) — C\ y\{x) = C 2 y 2 (x). 

If now A\ = we take ki = k 2 = [x] in (|1.14|) and get in a straightforward way: 
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1*1 „ m q 

°l,3 ^ __r —x \ "* °2. 



3=0 3=0 



< const x T 1//2 e 



Since by Lemma l3~Tl the same inequality is true with y(x) replaced by C\y\(x) + 
CiViix) it follows that \y(x) — Ciyi(x) — C 2 y2(x)\ decays faster tha.ny 2 (x) as x — > oo 
which is possible only if it vanishes identically. □ 

Proof of Lemma 11.11 

For this purpose we only need the recurrence (|3.17J1 . the asymptotic behavior of 
its solutions ()3.18j) and the condition (jl.K)j) . Assume we have a formal power series 
S '■= J2kLo SkX~ k that solves a canonical differential equation of the type considered. 
The recurrence (I3.17|) provides us in a straightforward way with a set of equations 
for aj,bj, j = 0,1,2,...: 

- sidj + s b j+2 = Tj for j = 0, 1, 2, ... (3.46) 

where Tj only depend on the series S and on the coefficients cti, h+2 with i < j. 

The equations ()3.46|) alone would not determine the aj, bi uniquely but the con- 
dition (jl.lUj) binds a; and bi together. The best way to see this is to divide (|3.17j) by 
s n , use (|3.18|) . ()1.10|) (and the assumption R ^ 0) to write an asymptotic expansion 
to all orders in - of the resulting equation and then equate the successive powers 
of A. 

n 

The upshot is the system: 

aodj — bj = Tj for j = 2, 3... (3.47) 

where Tj depend on the series S and on di,bi with i < j. But ao and a\ are 
determined from (|3.18|) . namely a = a, ai = — r, so that (|3.46|) . (|3.47p determine 
uniquely the a,, 6j. □ 



Proof of Lemma 11.11 

For simplicity we make a change of variables so that 

|a| = 1; x G 1R + 

The recurrence relation for the formal series solution is obtained in the usual 
way, by inserting the formal series in the differential equation, expanding everything 
out and identifying the powers of x: 
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(fc - 



fc=2 X p 

oo r oo r oo oo r /oo 

xi xi ^ xi " tJ \ ^ xi 

3=0 3=0 X 3=1 3=0 X \j=l X , 

oo j 

EjEUi-^ ( 3 - 48 ) 

fe=l x J 

where J consists of all the integer tuples (i, ki, k p ), ki > 1 and z > 0, with 
ki + k 2 + .. + k p + i = k. It follows, 

s/c = ~(k + r- l)s fc _i + V f Pli s kl ...s k (3.49) 
« Ji 

where r = + 2/ 2 ,o s i an d the index set Ji in the sum excludes the tuples with 
ki > k — 2 from J . We make the substitution (suggested by solving the linearized 
version of the recurrence) 

s k = a~ k T k+r r] k 

where T k+r = T(k + r) if 9ft(/c + r) > —1 and T k+r = 1 otherwise (again, to avoid 
unpleasant poles). For large k, 

V V 

rj k = + 2^ f P ,ir)k 1 ---Vk p a - — (3.50) 

Jl 1 k+r 

The existence of a formal solution of the differential equation follows from the 
existence of a solution to the recurrence relation ()3.50|) which is obvious. The first 
objective is then to show that the sequence {r] k } k converges; the only delicate step is 
to show that the sequence { \r] k \ } k is bounded. To this end we compare the recurrence 
(I3.5U|) with a suitable linear recurrence (|3.6Uj) below. 

Taking out of the sum the terms with ki > k — n for a conveniently large n we 
get, for some constants A\, .., A n , 



Vk = Vk-i + + •• + T^Vk-n + 2^ fprfki-Vkpa ^ (3.51) 

where the restriction J 2 in the sum now reads k-y + k 2 + .. + k p + i = k; ki < k — n. 
Taking the initial condition rji = \rji\ in the recurrence 
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Vk = Vk-i + -T^Vk-2 + ■■ + -T^Vk-n + 2_, f P ,iV kl -V kp r — IP = ^M) 

K J2 k+ p 

(3.52) 

we clearly get r] k > \i] k \. Now the solution of (J3.52J) with positive initial conditions 
is increasing in k so that for A > max{Ai, .., A n } 

Vk <( 1 + ^) + E UM^-v J^^ (3-53) 
Taking r] k = a k Y\\(l + A / j 2 ) / - 1/j 2 ) we get for some constant C 

*k< (1 - ±) + cY.U^k 1 -Vk p Tkl X' Tkp+P ( 3 - 54 ) 

Noting that card(J2)< P k p we get 



*k<(l- -jj) + ^ max |^+ 2 f Pli * kl ...* kp Tkl+ f^ kp+p I (3.55) 

It is easy to see by a straightforward calculation that one can choose k and n 
large enough so that 



-2* r fcl+p ...r fcp+p / k x \..k p \i\ y ( h\..k p \ 



CP fcP + V P ,i < „ ■ 7 .„ < „ r P \ „ (3.56) 



(A;i + .. + k p + i)\l \(fci + .. + fc p )! 



when ki < k — n (for the second inequality above note that the middle term is a 
decreasing function of its arguments). Using (jl.24j) we get 



° k ~ i 1 ~ ¥) ak ' 1 + ^T {(7fa -^ ((fc 1 + .'.' + ! A; p )!) 2} (3 ' 57) 

J 3 consists of all the tuples kx, .., k p with k\ + .. + k p < k. By construction the empty 
tuple contributes with const n h < const. It follows easily that 

a k < m f {^ i; a kl ...a kp ( ^f^^y) } (3-58) 
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from which, if G\ is large enough which we are allowed to assume since the solution 
is increasing in a\, it follows that a k is bounded by the solution of the recurrence: 



k\\..kp: 



1/2 



ak = ^ x{ak ^{(k 1 + '.'.+k p )l) } (3 - 59) 

For each k, ()3.59|) determines a tuple &*,..,&*, k\ + .. + k* < A; for which the 
maximum is realized. If we take a k = exp(0&) we get a linear recurrence: 

Consider the associated homogeneous recurrence 

F k = F K + ... + F k; - F 1 = l (3.61) 

Inductively we see that: 

F k < k (3.62) 
1 



We claim that 



(t>k<<PiF k --\n{F k \) (3.63) 

Indeed, the inequality above is true for k = 1 and assuming it holds for all k! < k 
we have (cf. also the comment following ()3.56j0 

in view of f)3.62|) which gives ()3.63|) . 

But note that the function g(x) = Ax — B ln(x!) : (l,oo) i— > H is bounded 
from above so that the solution of (J3.5(J|) is bounded. 

Now, taking n large enough one can make the nonlinear term in (J3.51j) smaller 
than const k ~ m (cf. also ()3.56^ with m as large needed. The rest of the proof is 
obvious. □ 

We return to the proof of Theorem 11.21 

Proof. We start by showing the existence of a special solution corresponding to 
C = 0. We look for solutions of the equation (jl.23j) in the form 

n Si 

y(x) = Y n (x) + V(n;x); Y n := £ -|; (xe[n,n + l]) (3.64) 

k=i x 
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where n = [x]. The differential equation for V reads: 



d p d p 

V = —F(Y n , x)V+J2 T^ F ( F - x ) VP + T - ( 3 - 65 ) 

dy ^2 d v p 



where 



T n = F(Y n , x) — Y' n (x) 
to which we add the condition of continuity of y: 



V(n,n + ) = V(n-\,n-)-?± (3.66) 



Lemma 3.2 There exists a solution of the problem hS. 6<ty - lS. 66]) with the property 

\V(n, x)\ < const n p (x G [n, n + 1]) 
Again we need first some estimates on the power series solution. 

Proposition 3.4 For x E [n,n+l] 

T n = (R + 0&)a-^ T -^ 



n x" 



Proof 

For p G {0, .., P} we write 



f p( x ) = £ H + with < C F n n n\ (x e[n,n + l}) (3.67) 

fc=0 x x 

Then T n is given by 



p=0 \k=0 x x / 

1 n k Sk P Rpn 

Zq fp,i s ki---Sk p + 22 ~XfT + z£n Y n 

Q X Jq fe=l X i=l X 

where the index set Jq contains tuples (z, k±, .., G {0, .., n} with z + k\ + .. + fc p = 
Q. The recurrence relation for Sk is such that all the terms with Q < n have to 
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compensate each other whereas for Q = n + 1 they add up to as n+ ix n 1 . We thus 
get for T n 



rp s n+l 



< E iEIA.ilKI-KI + E^ra (3-69) 



)>n+2 x J Q i=l x 



The last term on the rhs of ()3.69j) can be estimated by const T n+p+ i x n 2 . For 
the first term we write 

E + E ) ^J2\UM-K\ (3.70) 

yn+2<Q<n+s n+s<Q<(P+l)n/ Jq 

For all n, |s n | and |/ p , n | are bounded by const T n+P . Note that the T function 
is log-convex so that the maximum of a product of the form T(xi)..T(x m ) over a 
convex domain, is reached on the boundary of the domain. Now, if n is large enough 
compared to s + p, 

p+s+l 

T kl+p ..T kp+p < const P n ~ J E F fci 

ki + ..+k p =n+s 3=0 ki+..+k p -i=s+j 

ki<n ki<n 



n + p + S 



const ^ £ ^n~ V < const (P,s)r n+p (3.71) 

It follows that the leftmost sum in ()3.7()j) is less than const T n+p x~ n ~ 2 for x G 
[n, n + 1]. In the second sum we write 

Q = qn + m 

with q < P,m < n. The log-convexity argument together with the restriction on 
the indices i, kj give 

|/p,*||sjfcil — |sfc„| < const K l i\T kl+p ..T kp+p < 

const {T n+p ) q T m+p (3.72) 

and we crudely bound the second sum in ()3.7()|) by the number of terms times the 
maximal term. An elementary application of Stirling's formula gives 

constn p + 2 max ( BtgEgtg l < (3.73) 

n +s<Q<nP 1 n qn+m I V > 
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const n p+2 exp ( max {(m + p — -) ln(m) + (qp — ^ — m) ln(n) — m — gn}^ 

\n+s<Q<nP L 2 2 J / 

The expression to be maximized is decreasing in q and is convex in m. It is not 
difficult to see that for n large enough the maximum is reached at q — l,m — s. 
For large n (J3.73j) is bounded by const r f- s - 1 / 2 + p + 2 e _n . Choosing s > P + 4 we 
make the second sum in (J3.7U)) less than Y n+p n~ n ~ 2 . The proof is completed by 
combining these inequalities with (J3.69|) □ 

We take now n large enough and first find a suitable solution of the linearized 
version of the equation ()3.65|) on each interval [n,n + 1]. Consider the differential 
equations: 

H f = F y (Y n ,x)H + T n (3.74) 

with the initial conditions at x = n chosen so that y(x) as defined in ()3.64|) is 
continuous namely, 

H(n,n+) = H(n-l,nJ)-^ (3.75) 
under the same hypothesis on F as before. 



Proposition 3.5 There is a solution to the problem \3. 74 ), \S. 75] ) such that 



\H(n,x)\ < const n+p [x G [n, n + 1]) (3.76) 
x n 

The proof is the one- dimensional projection of the proof of lemma l3~D □. 
It also follows that for x £ [n, n + 1] 

H(n, x) < const x p ~ 1/2 e~ x (3.77) 

We now turn to the proof of Lemma 13.21 We shall look for solutions of the 
equation ()3.(i5j) in the form V = H + 5 where H is the function provided by the 
Proposition 13.51 Then S is a smooth function and satisfies the differential equation 

d p d p 

5' = —F(Y n , x)6 + J2 -^~ p F(Y n , x) (H + 6)? (3.78) 

We are looking for a solution of ()3.78|) of the order 0(H 2 ); such a solution is 
constructed as a fixed point of a contractive mapping. 
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Take n large enough, (3 > 2p — 1 and let S n ^ be the space of continuous functions 
on (n, oo) with the norm 

11/11 :=sup 

t>n 

and the operator 

r°° d p d p 

(Cf)(x):= j x -F(Y n ,t)f(t) + ^—F(Y n ,t)(H(t) + f(t)ydt (3.80) 

In view of (13.77)1 for large n, H 2 e «S n ,/3- Actually we can make the norm of H 2 
very small by choosing n large enough. L maps S n $ into itself and for a given e 
there is n e large enough so that C maps B e , the ball of radius e, into itself. Actually, 
£ is a contraction for large n. Indeed 

Cf= af(t)dt + 0(f/x,f 2 ,Hf,H 2 ) (3.81) 

Jx 

where O is with respect to the norm ()3.79|) . But if /, g e B t and e' > it is 
straightforward to see that \\Cf — Cg\\ < ^f-\\f — g\\ if n is large in which case 
the equation ()3.78|) has a solution in <S n>l g. It follows that the equation ()3.65|) has a 
solution with the required properties □. 

Lemma 13.21 together with Proposition 12.11 prove the existence part of Theorem 
11.21 for C = 0. Given this, for C / the existence result follows easily, since 
for any C there exists a unique solution of the differential equation (jl.23|) of the 
form y(x) = V(x) + e(x) with V(x) the solution given in Lemma 13.21 and e(x) ~ 
Ce~ ax (x r+l + 0(x r )) (\x\ — ► oo). This is a standard result (see [H], P) ; in our 
case e(x) could be also directly obtained using the contractive mapping arguments 
above. 

The uniqueness proof is very similar to the one in the linear case since no two 
solutions can differ by less than const e~ ax x p+1 as x —>■ oo without being equal to 
each other. Finally, the information that a solution y(x) is decaying at infinity 
is enough to guarantee that y is asymptotic to the series f)1.2fij) . The difference 
y(x) — V(x) has the asymptotic behavior Ce~ ax (x r+1 + 0(x r )) for some C. Let yc 
be the special solution satisfying (|1.30|) . Since y(x) — V(x) ~ Ce~ ax (x r+1 + 0(x r )) 
as well, it follows that y(x) — yc{x) = o(e~ ax x r+1 ) which means y(x) = yc(x) ( 9J; 
a simple direct proof uses contraction mapping ()3.80|) ) . The proof of Theorem 11.21 
is complete. □ 



f(t) 



(3.79) 
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Discussion of the examples. 
i) Proof of Eq. (jl.7|) . We are looking for a constant C, dependent on the ray in the 
complex plane, such that 



E n (x) := £i(x) - e x ^ ^ ^ = C + o(l) (\x\ -> oo, arg(x) = 0) (3.82) 

fc=i x 

By Theorem 1.1 and Proposition 12.11 there exists such a constant. We might as 
well compute it for \x\ — » oo along a subsequence, say x n = ne* 9 . 
• Take first 9 > (6* < is similar). It is easy to show that 



EJx) = m + n!n n dt 
Jc 



Ic t n+1 

where C is a contour joining — oo to e above the real axis; we choose it to be 
the stationary phase (or, which is the same, the steepest ascent) contour for the 
integrand: ^sint — in + 1) ln(£)) = const. 

The leading behavior of the integral for large n is due to the contribution of a 
region, near the right end point of C, where z := t — e td is of the order 0(l/n). 



f dt^-r = (1 + o(l))e ne ' B ~^ ie [° dze n ^- e " B > 

JC t n+1 K K " J-oc 



so that, for 9 > 0, 




E n - vrz = (1 + (1))^_ —e n ^ = o(l) (3.83) 



For 6* = we have: 



dt— T = -n!n" n e n S / 



-dt 



iyY+l 

where the last expression is obtained after pushing the contour parallel to the imag- 

1 _l^ 

2' 3' 



inary axis and taking t = 1 + iy. Now, choosing (3 G (— h, —I) we get 



-my 



iy) 



n+l 



dt 



-my 



iy) n 



+i 



dt + 0(exp(--n 2/3+1 )) 



e"^ (l+iy- '-ny^ dy + o^ 1 ) = ^- + JL + opi8i) 
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so that, after taking the imaginary part, 

E n {n) = -\J— (l + o(l)) (3.85) 
3 V n 

□ 

For the exponential integral Ei(x), it is easy to evaluate the optimal weight 
from ()3.85|) and its differential equation. The arguments below can be easily made 
rigorous but we will not insist on that, since we are only aiming at an illustration. 
Let 

n-2 u\ 

R(x) = e- x Ei(x) - E ^TT ( 3 - 86 ) 

fc=0 x 

Then, R(x) satisfies the differential equation 

(n-l)\ 



R' + R 

x" 

Taking R(x) = ^ n ~^' g(x) we get for g the equation g' + (1 — nx~ 1 )g = 1. The 
relevant region is x = n + 0(^/n) so it is convenient to change variables further to 
x = n + Sy/n and g(x) = y^n^(s). We get 



with the initial condition, coming from (J3.85j) 0(0) = 0(?t, _1 / 2 ). The solution is 
given by 



e 2 



[ S e^dt + 0(n- 1/2 ) 
Jo 

The maximum value of \(f>(s)\ is therefore equal to a*. 

ii) Derivation of (jl.21j) . To make the calculation of the asymptotic series easier 
we take further g(s) = s~ 1 ^ 6 h(s) (to make r\ = 0) and finally s = 3t/4 (to have 
Ai = 0; A2 = —1). The resulting equation is 

h" + h'+ -^h = 
36t 2 

having the following general formal solution: 
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00 h< „ + ^ f-)^. 



The recurrence for 



(j + i)Vi = (j+g)(j + ~)^ 



with the choice h = 1, gives 



= r(n + |)r(n + |) = _i_ / IN j 

27rr(n + 1) 2vr ; V nV 

so that with this choice we get R\ — Ri — ^ and /& = f|j 3 /i^. The contents 
of (jl.21|) is given by Theorem ll.il after making the appropriate substitutions and 
change of variable back into eq. ()1.14j) . 

4 Convergence of asymptotic series and a decom- 
position property 

In this section we give an example of of non-generic convergence of the 

asymptotic series and discuss its relevance for asymptotics beyond all orders. The 
Proposition below could be easily generalized to a larger class of differential equa- 
tions but we now look, for the sake of simplicity, at equations of the form 

ijj' + ijj = R{x) (4.1) 

where R(x) is a rational function. We consider ()4.1|) for large x on a ray arg(x) = 
9 G (— f , §)■ If R(x) is a polynomial the equation has an explicit solution of the 
form Polynomial(x) +Ce~ x so the interesting case is 

m 

R(x)= J2 R k>j (x- rj r k (4.2) 
k,j=i 

Proposition 4.1 i) Given R(x) of the form fl^.ifl ) there exists a unique constant K 
such that the solution of the equation 

K 

$ + tfj = R(x) (4.3) 

x 
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is holomorphic in a neighborhood of infinity. Consequently, any solution of the 
equation \4-l\j has a decomposition: 

ip(x) = K e~ x £i{x) + h(-)+C e~ x (4.4) 



where the constants K, C as well as the analytic function H(-) are uniquely deter- 
mined. 

ii) Let 9 = 0. A solution of \4-l\) has (in the sense of Theorem \l.ty) the asymptotic 
representation 

^ x )^J2- k +Ce~ x (4.5) 

3=0 X 

Moreover the constant C is the same as in \4-l\j and the power series is (factorially) 
divergent unless K = 0. 

Comments. On the one hand the proposition above indicates in what sense diver- 
gence of the asymptotic series is generic. On the other hand the decomposition (|4.4|) 
suggests another point of view on the problem of the terms beyond all orders for 
the equation (J4.1j) . Since in any reasonable definition of asymptotic representations, 
a function which is analytic at infinity should be represented by its own (conver- 
gent) asymptotic series, once the terms beyond all orders for a particular function 
(the exponential intergal) are defined, they can be determined unambiguously for 
the solutions of (J4.1)) with any rational inhomogeneity R. Part ii) shows that the 
results obtained in this way are consistent with those obtained through asymptotic 
estimates. 

Proof 

We have 

R( X ) = / e~ tx £ R k ^ t J ——dt (4.6) 

kj=i 



Let 



m -j 



Then the function: 
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is entire so that the function 

roo _ 

z / e- l R(zt)dt (4.9) 



o 

is analytic in z for \z\ < -. Indeed, the integrand is analytic in z and, \R(zt)\ < 

z m t m e' 2 ' p * so that the integral is uniformly convergent for \z p\ < 1. Furthermore, 
the function: 

1 r°° ~ f t \ 

- e- t R(-)dt + K£i(x) (4.10) 
x Jo \xj 

is a solution of the equation 

K 

$ + tfj = R( x ) (4.11) 

x 

as it can be easily checked. The rest of the proof of i) is immediate. 
Let 

OO 1 

Since clearly — > as n — * oo the proof of ii) is an easy application of Proposition 
Oand of formula (fT77j) .n 
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